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ABSTRACT 


The  resource  leveling  problem  for  a  construction 
system  producing  a  stream  of  output  units  Is  con¬ 
sidered.  The  system  Is  modeled  using  a  critical- 
path-analysis  activity  network,  from  which  an 
extended  network  is  developed  for  an  integrated 
planning  effort  of  all  output  units.  Activity 
intensity  variables  are  defined  which  measure 
activity  demand  rates  for  resources  and  consequent 
activity  durations  for  the  production  of  each  out¬ 
put  unit.  A  heuristic  approach  consisting  of  an 
iterative  non-linear  programming  procedure  is  pre¬ 
sented  which  computes  activity  durations  (intensi¬ 
ties)  for  the  minimization  of  resource  capacity 
costs  subject  to  meeting  construction  due  dates. 
The  application  to  a  major  ship  overhaul  is  de¬ 
scribed,  in  which  the  procedure  was  used  to  level 
workloads  of  the  various  labor- trade  shops. 
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Multiple  Resource  Leveling  in  Construction  Systems  Through 
Variation  of  Activity  Intensities 

Robert  C.  Leachman 
University  of  California.  Berkeley 

1.  Introduction 

In  construction  systems  such  as  shipbuilding  or  aircraft  assembly,  one  finds  the  construc¬ 
tion  of  a  stream  of  large  output  units  whose  resource  requirements  interface  through  time.  Typ¬ 
ically,  there  are  various  progress  and  completion  due  dates  on  each  output  unit,  and  the  prob¬ 
lem  arises  of  minimizing  peak  demands  for  various  resources  subject  to  meeting  these  due 
dates. 

In  the  usual  formulation  of  the  resource  leveling  problem  for  a  construction  project,  a 
critical  path  analysis  of  the  project  is  made,  in  which  component  activities  have  fixed  durations 
and  fixed  rates  of  resource  utilization.  The  problem  considered  is  to  schedule  the  activities  so  as 
to  minimize  resource  loading  peaks.  The  difficult,  combinatoric  nature  of  the  problem  has 
resulted  in  an  emphasis  on  heuristic  techniques  [21  [5]  (6!  [10], 

However,  in  many  construction  systems,  activities  are  "rubber  band"  in  nature;  that  is, 
there  is  a  feasible  range  of  intensity  of  resource  assignments  to  each  activity,  resulting  in  a 
range  of  possible  activity  durations.  In  this  paper,  the  leveling  of  the  various  resource  histories 
through  variation  in  the  intensity  of  activity  resource  loading  is  investigated. 

2.  The  Model 

2.1.  Extended  Network  For  Output  Streams 

We  consider  a  system  with  component  activities  i—l . /V  utilizing  K  exogenous 

resources  to  produce  M  output  units.  We  suppose  the  construction  of  an  individual  output  unit 

may  be  represented  by  a  CPM  node  scheme  activity  network,  in  which  the  nodes  represent  the 

activities,  and  the  arrows  indicate  work  flow  dependencies.*  The  network  is  to  be  sufficiently 

•  In  most  presentations  of  CPM/PERT.  (he  roles  of  nodes  and  arcs  are  reversed  from  the  roles  used  here 
Both  formats  are  used  in  practice.  See  (6|. 
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general  so  as  (o  model  any  of  the  output  units,  where  work  requirements  may  vary  from  output 
unit  to  output  unit. 

Sherrard  and  Mehlick  [?]  have  used  such  an  approach  to  perform  PERT  simulations  of  a 
series  of  construction  projects.  In  their  paper,  the  projects  are  performed  according  to  a  basic 
activity  network,  but  activity  durations  vary  from  project  to  project  due  to  varying  work  require¬ 
ments  and  learning  effects.  They  demonstrate  the  effects  on  completion  times  of  individual  pro¬ 
jects. 

However,  to  dynamically  model  resource  loading,  an  integrated  planning  effort  is  required, 
capturing  the  interaction  of  resource  demands  among  the  various  output  units.  We  consider  an 
extended  network  whose  nodes  are  denoted  by  im,  where  /  denotes  the  activity  and  m  denotes  the 
output  unit.  The  immediate  predecessors  of  node  im  in  the  extended  network  are  defined  as  fol¬ 
lows: 

(i)  <m_|,  except  when  m-l; 

(ii)  jm,  for  all  activities  j  which  are  immediate  predecessors  of  i  in  the  activity  network. 

Using  the  extended  network,  we  envision  the  construction  of  M  output  units  as  a  single  project. 

2.2.  Activity  Intensities 

The  dynamic  activity  analysis  of  production,  in  which  activity  output  rates  and  resource 
loads  are  measured  in  terms  of  intensity  functions  and  technical  coefficients,  is  due  to  Shephard 
et.  al.  [Si  [91.  These  authors  developed  the  analysis  for  continuous  flow  production  systems 
governed  by  linear  inequalities.  We  adapt  it  here  for  the  extended  network  format. 

The  rate  of  application  of  resources  to  each  activity  i  working  on  each  output  unit  m  is 
modeled  linearly  in  terms  of  a  strictly  positive  intensity  variable  :,m  taken  with  technical 
coefficients  defined  as  follows: 

(a)  cim  -  fraction  of  output  unit  m  completed  by  activity  /  per  unit  intensity  of  >  working  on 
m,  per  unit  time,  i-l . N .  m-I . A/. 

(b>  aimk  -  amount  of  exogenous  resource  k  applied  to  activity  /  per  unit  intensity  of  >  working 
on  output  unit  m,  per  unit  time,  /■*! . N,  m-1 . M,  A»1 . K. 

These  technical  coefficients  vary  from  output  unit  to  output  unit,  reflecting  varying  work 
requirements  as  well  as  learning  curve  behavior. 
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The  intensity  assignment  z,m  to  node  im  of  the  extended  network  is  maintained  at  a  con¬ 
stant  value  between  start  and  finish  times  for  im,  but  the  particular  level  of  intensity  is  a  deci¬ 
sion  variable.  In  this  way.  the  time  assignment  to  each  node  i„  becomes  a  variable,  given  by 

— — - — ,  /—I . N.  m—l . M .  (l) 

Z,mCim 

For  each  node  im  of  the  extended  network,  there  is  a  maximum  intensity  assignment 
termed  an  intensity  bound.  These  bounds  reflect  available  workspace  and  other  conditions  pecu¬ 
liar  to  each  activity. 

An  initial  selection  of  node  intensities  allows  the  ordinary  CPM  scheduling  computations 
(see  |6))  of  early  (late)  start  and  finish  times  and  total  activity  slack,  using  the  node  durations 
(1).  In  particular,  operating  all  activities  at  their  intensity  bounds,  ie. 

1 . -V.  m-1 . M.  (2) 

yields  a  minimum  project  time  span.  We  shall  refer  to  such  a  project  policy  as  the  greedy  plan. 

The  intensity  :,m  is  said  to  be  critical  if  the  time  assignments  given  by  (1)  are  such  that 
node  im  is  on  a  critical  path  in  the  extended  network.  If.  on  the  other  hand,  node  im  has  slack, 
it  is  an  indication  that  activity  i  working  on  output  m  could  be  operated  at  a  lower  intensity, 
using  less  resources  per  unit  time,  without  slipping  the  project  schedule.  In  general,  the  total 
slack  of  a  given  slack  path  through  the  extended  network  can  be  allocated  among  the  activities 
on  that  path,  whereby  the  reduced  intensity  assignments  would  become  critical.  In  fact,  all 
nodes  can  be  made  critical  by  eliminating  all  slack  from  the  slack  subnetwork,  ie.  the  subnetwork 
of  slack  paths  of  the  extended  network. 

In  following  sections,  we  attack  the  resource  leveling  problem  using  a  technique  for  allo¬ 
cating  the  slack  of  the  slack  subnetwork  to  reduce  node  intensity  assignments. 

3.  R  mm  res  Leveling 

3.1.  Fermnlatlen 

We  suppose  the  various  due  dates  for  the  project  are  expressed  as  late  finish  times  for  a 
subset  of  the  nodes  of  the  extended  network.  It  is  assumed  that  all  terminating  nodes,  ie,  nodes 
with  no  followers,  have  such  dates.  For  simplicity  of  exposition,  it  is  also  assumed  that  all  initial 
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nodes,  ie,  nodes  with  no  predecessors,  have  early  start  time  0.  Denoting  the  set  of  initial  nodes 
by  /  and  the  set  of  terminating  nodes  by  D,  there  is  thus  a  maximum  time  length  allowable  for 
each  path 


through  the  extended  network  starting  with  a  node  aW|  c  /  and  ending  with  a  node  «  D. 
This  length  is  given  by  the  late  finish  time  of  /3m,,  which  we  denote  by  LFU 3,„2). 

For  the  resource  leveling  problem,  we  consider  the  non-siorable  resources  required  by  the 
construction  system;  these  are  the  resources  which  cannot  be  accumulated,  such  as  the  services 
per  unit  time  of  various  machines  and  labor.  It  is  assumed  that  such  resources  cannot  be  "hired 
and  fired"  with  the  variation  of  project  resource  loading  histories,  and  consequently  resources 
adequate  to  meet  peak  demands  must  be  maintained  for  the  life  of  the  project. 

Let  NS  <  K  be  the  number  of  non-storable  resources,  and  renumber  (if  necessary)  such 

resources  4-1 . NS.  Assuming  linear  capacity  costs,  let  pk  be  the  cost  of  maintaining  a 

unit  of  capacity  of  resource  4.  4-1 . NS,  for  the  life  of  the  project,  where  such  capacity  is 

required  to  meet  peak  demands. 

To  simplify  the  calculation  of  project  resource  loads,  we  consider  a  bar  graph  approxima¬ 
tion.  We  take  a  discrete  time  grid  t  —  0, 1 , 2 . T.  where  T  is  the  last  project  due  date,  and 

the  grid  is  fine  enough  so  that 

. N  m-1 . M.  (3) 

That  is,  each  node  requires  at  least  one  time  unit  duration.  For  a  given  set  of  node  durations 

/„„J,  /-I . N.  m-1 . M. 

one  calculates  the  project  loading  of  resource  4  during  [r,r+ 1)  for  an  early  start  schedule  to  be 

(4) 

i-l  m-l  C"» 


where 

jl  if  r-t,m  <  £S(tm)  ^  t, 
#(t.m,r)-|o  if  no( 


(5) 
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Here,  ES(i„ )  denotes  the  early  start  of  node  The  delta  functions  defined  by  (5)  are  seen  to 
indicate  whether  node  im  operates  during  (r,  r+1)  or  not. 

We  consider  the  problem  of  minimizing  costs  for  non-storable  resource  capacities,  subject 

to  meeting  project  due  dates.  Let  Xk.  Ar — 1 . VS,  be  variables  representing  the  peak 

demands  for  non*storab!e  resources  encountered  by  the  project.  The  problem  may  then  be  for* 
mulated  as  follows. 

vs 

Minimize  ^.Pk^k 

4-1 

subject  to 

*<l>  I  tm  <  LF(pmi). 

'm  ""  F«„0„2 

for  all  /3,„  c  D,  and  all  a,„  c  /.  and  for  all  paths  Pa  a  through  the  extended  network  starting 

«  I  Wj 

at  am|  and  ending  at  f3my 

R( 2)  £  £  Sd.m.r)  —  —  -  XK  <  0.  A-I . NS.  r-0. 1 . T- 1. 

im\ 

where  is  defined  by  (5). 

/?  (3)  >  ■“ — .  /-I . N.  m- 1 . A/. 

‘im  Cm 

/?  (4)  /,m.  AT*  >0.  /-I . N.  m- 1 . A/,  k- 1 . AS. 

Constraints  A (2),  for  which  one  must  relate  node  start  times  to  the  time  grid 

r  -  0, 1 . r-1,  are  a  source  of  much  difficulty.  The  literature  contains  efforts  to  handle 

constraints  of  this  kind  using  integer  variables,  making  the  problem  too  large  to  solve  for  all 
but  very  short  time  horizons  T.  A  convenient  survey  of  such  efforts  is  Herroelen  (31.  We  shall 
resort  to  a  heuristic  approach,  as  described  in  the  next  section. 
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3.2.  Peak  Pricing  Procedure 

The  procedure  presented  here  arises  from  the  observation  that,  in  many  projects,  activities 
in  the  very  early  and  very  late  stages  of  the  project,  while  requiring  time  to  perform,  are  not 
resource  intensive.  Such  activities  might  include  tests,  inspections,  system  checkouts,  etc. 
Thus  resource  demand  peaks  for  such  projects  are  somewhat  inflexible  as  to  time  of  occurrence, 
restricted  to  mid-project  ranges  of  the  time  grid. 

A  sequence  of  problems  are  considered  in  which  only  nodes  operating  in  time  regions  of 
peak  resource  demand  in  the  solution  at  hand  are  charged  for  resource  capacity.  Following 
optimization  of  each  problem,  inspection  is  made  to  see  if  peak  regions  have  shifted,  and,  if  so, 
to  redefine  the  set  of  nodes  to  be  charged  for  capacity  in  another  iteration.  The  procedure  is 
presented  in  algorithmic  form  as  follows. 

Step  0.  Compute  the  greedy  plan,  ie,  set  node  intensities  and  time  assignments  as  in  (2)  and 
(1),  respectively.  Compute  a  CPM  forward  pass  schedule  for  these  time  assignments,  and  if  all 
project  due  dates  are  not  met,  stop;  the  project  is  infeasible.  Otherwise,  go  to  Step  2. 

Step  1.  Compute  a  CPM  forward  pass  schedule  for  the  node  time  assignments  of  the  current 
solution. 

Step  2.  Compute  the  project  resource  demand  histories  for  the  current  solution  using  (4)  and 

(5).  From  these  histories,  for  each  resource  Ar  -  l . VS,  identify  time  regions  of  peak 

demand  to  be  'priced,*  whereby  activities  operating  in  such  regions  will  be  charged  for  capacity 
in  this  iteration.  The  region  specification  is  somewhat  arbitrary,  but  should  include  time  inter¬ 
vals  experiencing  peak  and  near-peak  demand,  and  exclude  intervals  of  low  resource  demand. 
Time  intervals  excluded  from  peak  regions  are  likely  to  experience  higher  levels  of  resource 
demand  in  the  next  solution.  See  Figure  1  for  an  example  of  region  pricing. 

Based  on  these  peak  region  definitions,  compute 

1  if  8 (i.rn.r)  -  I  for  some  t  such  that 
bii.m.k)  —  Ir,  r+1 )  belongs  to  the  peak  region  of  k, 

0  if  not, 

i-l . /V.  m-I . M,  Ar—1 . NS.  (6) 


If  the  pricing  scheme  (6)  has  not  changed  since  last  iteration,  stop;  otherwise,  solve  the  follow¬ 
ing  program. 


Minimize 


1 


VS  V  V 

ill  b(i.m.k)  pk 


@imk 


K  ~  1  /  —  l  »»■  I 


I  “ 

'  *•/ 


+  5.., 


subject  to 


A/(l)  I  s,w  <  -  I  (F,„,c, 

/  on  P  _  ,i  /  „  on  P  ,» 


for  all  /3,„  «  A  and  all  a,„  t  /.  and  for  all  paths  Pa  0  through  the  extended  network  starting 

*  I  w|  mi 

at  a„t|  and  ending  at  /3m,. 


M (2)  sim  >  0.  /-I . N.  m-1 . M. 


Step  3.  Set  t,„  —  s,m  +  (zmc„„)  /—I . ,V,  w— 1 . U,  and  go  to  Step  1. 

The  procedure  is  seen  to  terminate  when  peak  regions  stabilize.  We  shall  refer  to  the 
final  set  of  node  time  assignments  generated  by  the  procedure  as  the  resource  leveled  plan. 

In  each  iteration  a  nonlinear  objective  is  minimized  subject  to  linear  constraints.  The  pro¬ 
gram  is  formulated  in  terms  of  slack  allocation  variables 

a  /-I . n.  m-i . m. 


so  that  it  may  be  conveniently  restricted  to  those  nodes  belonging  to  the  slack  subnetwork  for 
the  greedy  plan  calculated  in  Step  0.  Constraints  Mil)  would  then  effectively  consider  only 
paths  in  the  slack  subnetwork. 

There  is  no  particular  justification  for  starting  the  procedure  with  an  early  start  schedule  of 
the  greedy  plan;  a  late  start  schedule  could  also  be  used.  At  the  end  of  each  iteration  of  the 
procedure,  however,  all  nodes  are  made  critical.  This  can  be  seen  by  noting  that  any  slack  on  a 
path  from  an  initial  node  to  a  terminating  node  could  be  allocated  to  those  activities  operating 
in  a  region  of  peak  resource  demand,  thereby  decreasing  the  objective.  Thus  the  procedure 
looks  only  among  project  plans  in  which  all  intensities  are  critical. 

It  is  to  be  emphasized  that  the  procedure  is  a  heuristic,  and  it  is  not  difficult  to  construct  a 
simple  example  with  very  flexible  peak  region  locations  which  will  cause  the  procedure  to  cycle 
two  solutions.  However,  for  such  a  case,  the  pricing  scheme  (6)  may  be  readily  adapted  to  allow 
the  immediate  escape  from  such  a  cycle  by  simply  pricing  peak  regions  from  both  solutions. 
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3.3.  Applications  and  Computational  Experience 

The  Peak  Pricing  Procedure  was  applied  to  a  network  consisting  of  1 5 1  activities  modeling 
a  ship  overhaul.  Since  only  a  single  overhaul  was  being  modeled,  the  activity  network  served  as 
the  extended  network.  Exogenous  inputs  to  the  activities  consisted  of  labor  hours  of  twelve 
different  trade  shops.  Unit  capacity  costs  were  developed  from  average  wage  and  overhead 
rates  for  each  shop. 

Computation  of  the  greedy  plan  established  a  minimum  project  time  span  of  370  working 
days.  However,  such  a  plan  was  characterized  by  serious  demand  peaks  for  various  trade  shops. 
The  problem  of  reducing  the  costs  of  manpower  capacity  without  slipping  the  project  span  was 
addressed. 

The  slack  subnetwork  under  consideration  was  reduced  to  that  consisting  of  activities 
whose  start  times  and  slack  possibly  allowed  them  to  operate  in  the  middle  time  range  of  the 
project  where  resource  peaks  could  occur.  The  reduced  network  for  application  of  the  Peak 
Pricing  Procedure  then  consisted  of  33  activities  contained  in  48  slack  paths. 

The  computer  code  used  in  the  Peak  Pricing  Procedure  applied  the  Frank- Wolfe  lineariza¬ 
tion  algorithm,  which  computed  a  solution  within  a  1%  tolerance  of  the  optimum  for  each  sub¬ 
problem.  (See  (ll  for  a  description  of  the  algorithm.)  Peak  resource  regions  were  defined  to  be 
all  work  days  where  resource  demand  was  within  4%  of  peak  demand. 

After  four  iterations,  peak  regions  stabilized,  so  that  no  further  improvement  was  possible 
for  the  Procedure.  Each  iteration  had  required  about  10-15  seconds  effective  time  on  a  CDC 
6400  computer.  Figure  2  compares  work  loads  of  the  mechanical  shop  for  the  greedy  produc¬ 
tion  plan  and  the  resource-leveled  production  plan.  As  can  be  seen,  considerable  improvement 
in  shop  loading  has  been  made.  Overall,  a  38%  improvement  in  the  objective  function  was 
made  between  the  initial  greedy  plan  and  the  last  iteration  of  the  pricing  procedure. 

The  Peak  Pricing  Procedure  was  also  applied  to  a  two-output  unit  extended  network  of  42 
nodes,  developed  from  a  21  activity  subnetwork  of  the  aforementioned  ship  overhaul  network. 
Figure  3  compares  workloads  of  the  shipboard  mechanical  shop  for  the  late  start  schedules  of 
the  greedy  plan  and  the  resource  leveled  plan.  In  this  application,  a  similar  number  of  itera¬ 
tions,  computer  run  times,  and  objective  function  improvement  were  experienced  as  in  the  pre¬ 
vious  application.  It  appears  that  considerably  larger  problems  could  be  handled  by  the  pro¬ 
cedure. 
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3.4.  Intensity  Lower  Bounds 

Considerable  replanning  of  activities  to  lower  intensity  levels  was  made  by  the  pricing  pro¬ 
cedure  to  obtain  the  load  leveling  discussed  in  the  applications.  In  many  systems,  there  may  be 
practical  lower  limits  on  activity  intensity,  ie.  intensity  lower  bounds 

r,*,,  /— 1 . N.  m— I . M. 

An  important  special  case  is  an  activity  which  can  only  be  performed  "one  way,"  ie.  upper  and 
lower  intensity  bounds  are  equal. 

Such  limitations  can  be  incorporated  into  the  pricing  procedure  in  the  form  of  upper 
bounds  on  node  time  assignments.  The  constraints  A/(  I )— A/<2)  need  only  be  supplemented  by 
additional  constraints  as  follows. 

,1/(3)  s,„  <  m>~'  -  (:imc„ /-I . ,V.  m-I . M. 

In  this  situation,  all  node  durations  might  not  be  critical  in  the  solutions  generated  by  the 
Peak  Pricing  Procedure,  as  the  node  duration  upper  bounds  may  be  such  that  the  absorption  of 
all  slack  in  the  network  is  impossible.  If  such  is  the  case,  different  schedules  of  the  resource 
leveled  plan  may  produce  different  resource  loading  histories,  and  so  a  two-phased  procedure  is 
suggested.  First,  application  of  the  Peak  Pricing  Procedure  is  made  to  determine  the  resource 
leveled  plan.  Second,  application  of  an  heuristic  scheduling  algorithm  such  as  developed  in 
references  121.  (51  or  (101  could  be  made  to  search  for  an  improved  rescheduling  of  the  early 
start  schedule  generated  by  the  pricing  procedure. 

3.5.  Activity  Load  Leveling  Program 

In  the  special  case  that  no  two  activities  share  the  same  resource,  an  alternative  non-linear 
programming  formulation  allows  optimal  solutions  to  be  obtained  for  the  resource  leveling 

problem,  which  in  this  case  amounts  to  an  activity  load  leveling  problem.  Let  NS  (/)./-! . V. 

denote  the  number  of  non-storable  resources  required  by  activity  i.  Then  the  program  is  to 

Minimize  £  £  ft  7- 

,-l  4-1 

subject  to 

S(I)  I  s,m  <  LFWm2)  -  I  (r,mc,w)-'. 

•»*  im  P..  rt 
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for  ail  <  0.  and  all  aW|  <  /.  and  for  all  paths  ^  through  the  extended  network  starting 
at  a,M|  and  ending  at  &m;. 

S( 2)  —  B,k  -  s,m  (:,„c,mr'.  /-I . N.  m-1 

C,m 

S(3)  s,„.  Blk  >  0.  i-l . N.  m-1 . M.  k- 1 

Optimal  node  time  assignments  are  then  obtained  by  setting 

/,*  -  S,„  +  /-I . N.  m-1 . M.  (7) 

Here,  each  variable  B,k  is  the  inverse  of  the  peak  demand  by  activity  /  for  its  resource 

k.  Ar— 1 . i—l . Constraints  S(2)  define  these  variables.  The  problem  has 

been  formulated  with  these  variables  so  as  to  obtain  a  constraint  set  of  linear  inequalities,  mak¬ 
ing  for  ease  of  solution.  Other  variables  and  constants  are  defined  as  in  the  Peak  Pricing  Pro¬ 
cedure.  Lower  intensity  bound  constraints  A/(3)  may  also  be  introduced  into  this  problem. 

4.  Conclusion 

It  has  been  shown  how  critical  path  analysis  may  be  generalized  into  a  production  planning 
and  scheduling  tool  for  construction  systems.  The  Peak  Pricing  Procedure  (or  Activity  Load 
Leveling  Program)  applied  to  an  extended  network  can  be  used  to  plan  resource  capacities  and 
schedule  activity  operations  subject  to  given  production  due  dates. 

The  extended  network  model  of  a  system  producing  a  stream  of  output  units  constitutes  a 
discrete  transfer  model  of  production.  That  is.  completed  units  of  work  by  an  activity  which  are 
required  as  inputs  by  other  activities  are  transferred  unit  by  unit  on  an  event  basis.  The  consid¬ 
erable  detail  provided  in  modeling  work  in  process  in  this  fashion  may  prove  to  be  intractable 
unless  these  units  are  'large,'  as  is  the  case  for  the  shipbuilding  example  which  motivated  this 
model.  In  models  of  other  kinds  of  production,  a  considerable  number  of  output  units  are  pro¬ 
ducible  per  unit  time,  and  transfers  of  work  in  process  may  reasonably  be  viewed  as  continuous 
flows  governed  by  linear  inequalities,  as  in  references  [8]  and  (91.  A  forthcoming  paper  by  the 
author  will  discuss  production  planning  with  continuous  flow  models. 
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